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Abstract 

This paper is devoted to the stabilization problem for nonlinear 
driftless control systems by means of a time-varying feedback control. 
It is assumed that the vector fields of the system together with their 
first order Lie brackets span the whole tangent space at the equilib¬ 
rium. A family of trigonometric open-loop controls is constructed to 
approximate the gradient flow associated with a Lyapunov function. 
These controls are applied for the derivation of a time-varying feed¬ 
back law under the sampling strategy. By using Lyapunov’s direct 
method, we prove that the controller proposed ensures exponential 
stability of the equilibrium. As an example, this control design pro¬ 
cedure is applied to stabilize the Brockett integrator. 


1 Introduction 

Consider a control system 

m 

X = Ef Uifi( x ) = f(x,u), (1) 

i=l 
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where x = (aq, x 2 , ..., x n )' G D C K n is the state and u = (tq, w 2 , u m )' G 
M m is the control. The domain D contains the trivial equilibrium point x — 0. 
We treat all vectors as columns and denote the transpose with a prime. The 
vector helds fi(x) are assumed to be mappings of class C 2 from D to M n . 

ft is a well-known fact due to R.W. Brockett |9j that system ([Tj) is not 
stabilizable by a smooth feedback law u = k(x) such that k( 0) = 0, provided 
that m < n and /i(0), / 2 (0), ...., f m ( 0) are linearly independent vectors. 
Note that Brockett’s condition remains necessary for the stabilizability in a 
class of discontinuous feedback laws provided that the solutions of the closed- 
loop system are defined in the sense of A.F. Filippov [2Tj. To overcome this 
obstruction, two main strategies can be used for the stabilization of general 
controllable systems. The first strategy is based on the use of a time-varying 
continuous feedback law u = k(t,x ) to stabilize the origin of a small-time 
locally controllable system |12j . In the other strategy, the equilibrium of an 
asymptotically controllable system can be stabilized by means of a discontin¬ 
uous feedback law u = k(x), provided that the solutions (“7r-trajectories”) 
are defined in the sense of sampling mi. 

An approach for the practical stabilization of nonholomomic systems 
based on transverse functions is proposed by P. Morin and C. Samson [18] . 
A survey of feedback design techniques is presented in the book by J.- 
M. Coron [12]. Despite the rich literature in this area and to the best of 
our knowledge, there is no universal procedure available for the stabilizing 
control design for an arbitrary nonlinear system of form (JT]). 

The paper [TO] is devoted to the control design for a kinematic cart model 
with two inputs. A coordinate transformation from the three-dimensional 
state space to a two-dimensional manifold (parameterized by the arc length 
and the orientation error) plays a crucial role in the analysis. Based on 
this representation, a discontinuous feedback law is proposed such that any 
solution of the closed-loop system exponentially converges to an equilibrium 
point. The orientation angle is defined modulo 2n in such equilibria. 

Applications of sinusoidal controls to the steering problem for systems of 
form (JT]) are considered in the paper [TU] . A combination of constant controls 
and sinusoids at integrally related frequencies is used to steer the first-order 
canonical system to an arbitrary configuration. Some modifications of this 
algorithm are presented for chained systems. An overview of algorithms for 
the motion planning of nonholonomic systems is presented in the book [15] . 

In the paper [I], the controllability and trajectory tracking problems are 
considered for a kinematic car model with nonholonomic constraints. A result 


2 


on the solvability of the motion planning problem is established for such a 
model by using trigonometric controls. The error dynamics in a neighborhood 
of the reference trajectory is studied to solve the tracking problem. It is 
shown that the error dynamics is stabilizable by using a quadratic Lyapunov 
function. The controller design scheme proposed is illustrated by examples of 
a state-to-state control and tracking a circle with time scheduling at selected 
points. 

The stabilization problem for a nonholonomic system in power form with 
bounded inputs is considered in the paper |3j. The receding-horizon principle 
is used to solve an open-loop optimization problem and to derive a sampling 
control. It is proved that the family of controls obtained can be used to 
stabilize the destination state in finite time with any chosen precision. The 
numerical implementation of this algorithm is shown for a five dimensional 
system. 

The paper [2Tj is devoted to the stabilization problem of nonholonomic 
systems about a feasible trajectory, instead of a point. For such kind of 
problem, a time-varying feedback law is obtained by using the linearization 
around a feasible trajectory. The Heisenberg system and a mobile robot 
model are considered as examples for stabilizing a straight line trajectory in 
the three-dimensional space. This approach is shown to be applicable for the 
trajectory stabilization of a front wheel drive car. 

Assume that m < n and that /i(x), f 2 (x),..., f m (x) together with a fixed 
set of the first order Lie brackets span the whole tangent space for system (JT() , 
i.e. 

span {fi(x), [fj,fi]{x) \i = 1,2, ...,m, (j, l) G S} = R n , (2) 

for each x G D, where S C {1, 2,..., m} 2 , 




dfi(x ) 
dx 


fj(x) 


dfj(x ) 
dx 


fii?) 


and d-bL) j s the Jacobi matrix. Without loss of generality, we assume that 
each pair (J, l) e S is ordered with j < l. 

Following the idea of [221 23], we introduce an extended system for (pQ): 


x = 5>/ 4 (*) + ^ UjilfjJiKx ) = f(x,u) (3) 

i=i (j,i)es 


with the control u = (iq, u 2 -,u m , u ji)'^ n eS - Because of the rank condi¬ 
tion (]2D, every smooth curve is a trajectory of system ( 0 ). As subspaces 
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spanned by the Lie brackets of vector fields f 3 (x) play a crucial role in the 
dynamics study of system (JTD , we note that harmonic inputs naturally ap¬ 
pear as optimal controls implementing the motion along a Lie bracket [8], [13j. 
A result on the convergence of solutions of system fl| to a solution of ([3]) 
is established by H.J. Sussmann and W. Liu. ft is shown in the paper [22] 
that if a sequence of input functions {td (f)}°L 1 of class L 1 (0, r) satisfies cer¬ 
tain boundedness condition and converges to an extented input u(t) in the 
iterated integrals sense, then solutions x^(t) of system (CQ) with initial data 
x J '(0) = x° converge to a solution x°°(t) of system (13|) . uniformly with respect 
to £ € [0, r]. This result is stated for an extended system with higher order 
Lie brackets as well. The problem of approximating a given trajectory of 
the extended system by trajectories of system (P) is solved in the paper [17] 
by using an unbounded sequence of oscillating controls with unbounded fre¬ 
quencies. For a class of control systems with periodic solutions and small 
controls, an averaged control system is constructed in the paper [6]. It is 
proved there that solutions of the averaged system approximate all solutions 
of the oscillating system as the frequency of oscillations tends to infinity. 

In contrast to the above approach, we will use a time-varying feedback 
control u = u(t , x) with bounded frequencies to implement certain decreasing 
condition for a Lyapunov functon along the trajectories of system (pQ). The 
rank condition (J2]) implies that any positive definite function V(x) of class 
C 1 (M n ) may be taken as a control Lyapunov function for system ((2]), so its 
origin x = 0 is stabilizable by a smooth feedback law u = u(x), h(0) = 0. 
Suppose that such a feedback u(x) is given, then our goal is to construct a 
time-varying feedback law u = u(t,x ) for the original system ([Tj) in order 
to approximate the flow of the closed-loop system (J3]) in a suitable way. By 
exploiting this idea, we establish a result on the exponential stabilization in 
the sense of sampling controls and “7T £ -solutions”. 

We prove that, for systems satisfying the rank condition (|2|) , there exists a 
feedback u = u £ (t,x) such that any 7r e -solution x(t) together with u £ (t,x(t)) 
tend to zero exponentially, provided that £ > 0 is small enough (Theorems [Tj 
and [2] in Section [2]) . The proof of this result, given in Section [T] is based on 
Lyapunov’s direct method and the representation of solutions by means of 
the Volterra series described in Section [3] The construction of a stabilizing 
control u = u £ (t,x ) is carried out explicitly in Section [5] for the Brockett 
integrator. We show that such a feedback ensures exponential stability of 
the equilibrium. 
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2 Stabilization with sampling controls 

For a given e > 0, we denote by 7r e the partition of [0, +oo) into intervals 
Ij = \pji tj+ i), tj = J = 0; 1) 2,... . 


The following dehnition extends the notion of “7r-trajectories”, introduced 
in HU, for the case of a time-varying feedback law. 

Definition. Assume given a feedback u = h(t, x), h : [0, +oo) x D — y 
e > 0, andx 0 G IR n . A n £ -solution of system (JT]) corresponding to x° G D and 
h(t,x ) is an absolutely continuous function x(t) G D, defined fort G [0,+oo), 
which satisfies the initial condition x(0) = x° and the following differential 
equations 

x(t) = f(x(t), h(t, x(tj))), telj = [tj, t j+ 1 ), 
for each j = 0,1,2,... . 

In order to stabilize system ([I]), we will use a time-varying feedback con¬ 
trol of the form 



on each interval Ij of length £, where e; denotes the i-th unit vector in M m , 
and functions v(x) = (vi(x),v 2 (x), ...,v m (x))', au(x), ku{x) will be defined 
below. 

Note that there is no control Lyapunov function for the original system © 
due to Artstein’s theorem [2] and Brockett’s condition j9j. Even though a 
Lyapunov function may be constructed for system (CQ) in the sense of partial 
stability pa. such partial formulation is not sufficient to establish an expo¬ 
nential stability result. Because of the rank condition ()2]h any differentiable 
positive definite function V : D —> M is a control Lyapunov function for the 
extended system (J3J). Our main idea is to choose the feedback control (J3J in 
order to approximate the direction of —'VV(x) by trajectories of system (P), 
where VL(x) is the gradient of V(x). For this purpose, we fix x G D and 
£ > 0, and consider the following system of second order algebraic equations 


^Mx) + 

i= 1 



ViVj 


dfj(x) 

dx 


fi(x)+ 
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[fi,f j ](x) = -VV(x), (5) 


£ 

2tt 


E U E E 


i<j 


(q,i)es 


(qj)es 


"03 




with respect to variables tg, a q i, i G {1,2, ...,m}, (g, Z) G 5, assuming that 
the numbers k q i G Z \ {0} are chosen without resonances, i.e. 


IM ^ IM for all (q,l) G S', (j,r) G S, (q,l) ± (j,r). (6) 

Let us denote by B p ( 0) C M n the open ball of radius p centered at x — 0, 
and let B p ( 0) be its closure. In this paper, we use the standard Euclidean 

norms for all vectors and treat 9 as the Hessian matrix of the j-th 
component of f 7 fx). 

The basic result of this paper is as follows. 


Theorem 1 Let V(x) be a function of class C 2 (D) such that 


||VH(a;)|| 2 > ai V(x), V{x) > ^\\x\\ 2 , V(0) = 0, 


( 7 ) 


and let 

< L, \/x G D, i G {1, (8) 

with some positive constants a±, /3 i; and L. Assume that, for some p 0 > 0 
and e 0 > 0, algebraic system (0 admits a solution 


dfijx ) 
dx 


Vi = v/l(x), aji = a e jt (x), kji = k e jt (x), i G {1, (j,l) G S, 


defined for all x G S Po (0) C D and e G (0, £o]> such that condition (j6j) holds 
and 


lim 

£—>-0 


H^ £ (x)|l + ||a £ (x)|^ 2/3 

O Lip II 111 / o ° 

0<||x||<po Ill'll 


0. 


( 9 ) 


Then there exist positive numbers p < p 0 and £ < £q such that, for any 
£ G (0,e], there is a X = A(e) > 0: 


x° G B p ( 0) =>■ ||x(t)|| = 0(e xt ), \\u e (t, x(£))|| = 0(e Ai//3 ) as t —>■ +oo, 

( 10 ) 

for each 'K e -solution x(t) of system (JTD with the control u = u e (t,x) of 
form 0). 
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Property (fT2]) implies, in particular, that all 7r e -solutions x(t) of the 
closed-loop system (JTj) and (J3J) with initial data ||x°|| < p are defined for 
all t > 0. 

To ensure the local solvability of equations (J5J) in some A-neigborhood of 
the point x = 0, we use the following lemma. 


Lemma 1 Assume that the vector fields fi(x), / 2 (:r), ..., f m (x ) satisfy the 
rank condition (J2D in a domain D C M", 0 G D, l^l = n — m, and let 
V G C 2 (D) be a positive definite function. Then, for any small enough 
e > 0, there exists a A > 0 such that algebraic system (J5]) has a solution 


v £ {x) = (u?(ar), ...,v E m (x)y , a £ (x) = (a £ fix) W)e s)', k e (x ) = (k^fx)^ eS )', 
such that conditions dHJ) hold for each x G B& (0). The above solution satisfies 


estimates 


|v £ 0)|| < My||x||, ||a £ (a;)|| < M a 


\x 


x G 5a(0), 


( 11 ) 


where positive constants M v and M a do not depend on e. 


The proof of Lemma [T] is based on the degree theory and will be presented 
in Section [4] Lemma Q] allows us to formulate a local version of Theorem [T| 
as follows. 


Theorem 2 Assume that the vector fields fi(x), •••, f m (x) satisfy the 

rank condition (J2]) with l^l = n — m at x — 0. Then, for any positive definite 
quadratic form Vfx), there exist constants po A P > 0 and £q > e > 0 such 
that algebraic system (J5]) admits a solution 

Vi = v £ (x), ttji = djfix), kji = k £ fix), x G B P0 ( 0) C D,£G (0,e 0 ], 

i G {1,..., m}, (jj) G S , 
and, for any e G (0,e], there is a X = A(e) > 0: 

x° G B p ( 0) =>■ ||x(f)|| = 0(e~ xt ), \\u e (t, x(t))|| = 0(e _A ^ 3 ) as t —>• +oo, 

( 12 ) 

for each 7T e -solution x{t) of system (JT]) with the control u = u e (t,x) of 
form © . 

Proof. The assertion of Theorem [2] is a straightforward consequence of 
Theorem [H To ensure condition (J9j), we use inequalities (TITj) from Lemma [U 
□ 

The next section provides some technical results for the control design and 
stability analysis. Then the proof of Theorem [1] will be given in Section [3 
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3 Oscillating controls and representation of 
solutions 

Any solution x(t) of system ([1]) with initial data x(0) = x° and controls 
Ui = Ui(t), Ui G L°°[0,r] can be represented by means of the Volterra type 
series (cf. [71 [2D]): 

x{t) = x°+y 2 fi(x°) [ Uj(t)dt+\ '' - fj(x°) [ m(t)dt [ Uj(t)dt+ 

i=1 Jo Z i,j= i 0x Jo Jo 

+ \ fj]( x °) [ [ {uj(t)u i {s)-ui{t)u j (s)}dsdt +R(t). (13) 

2 i<j 

Here, and in the sequel, 1 stands for the Jacobian matrix of fj(x) evalu¬ 
ated at x = x°. The remainder R(r) of expansion ([17]) is estimated by using 
the following lemma. 

Lemma 2 Let D C W 1 be a convex domain, and let x(t) G D, 0 < t < r, be 
the solution of system (P corresponding to an initial value x(0) = x° G D 
and control u G C[0,r]. If the vector fields fi(x), ,/ 2 (x), f m (x) satisfy 
assumptions 


dfi(x ) 

< L, 

d 2 fij( x ) 

dx 


d 2 x 


<H, 


i — l,m, j — 1, n, 


(14) 


in D with some constants H, L > 0, then the remainder R(r) of the Volterra 
expansion A3 satisfies the following estimate: 


IlflWII < f; { eWT ~ \ V Ut + it + 4 } + 


HM r ,.lut 
4L 3 V 


2) 2 + 2Lf/r-l| 


M(L 2 + HMy/n) 
6 


f/ 3 r 3 + 0(f/ 4 r 4 ). 

(15) 


m 

M = max ||/j(a; 0 )||, U = max ^ K(t)|. 

1 <i<m 0 <t<r ' 


Here 


















The proof of Lemma [2] is given in Section [3J 

In order to use the control strategy (jj]), we consider a family of open-loop 
controls 


Ui(t) =Vj+ ^ a# 
U,i)es 




+ 5 a sin 


2lxk n ^ 



1 , 2 , m, 


(16) 

depending on parameters v = (ui, v 2 ,v m )' G M m , a = e 

fc = (kji)'(j ^ £S G (Z \ {0}) n_m , and £ > 0. Here 5^- is the Kronecker delta. 

By computing the integrals in (fT3j) for functions n* = Mj(t) given by (1T6T) 
and exploiting assumption ((HI), we get 


x(e) = x° + £ ^2 v iM x °) + 

i=1 


£ 

"2 


E 

*h=i 


ViVi 


dfj 


x 


dx 




+ < ^{ S jl( a qldiq'^ 2v i)- S il( a ql 5 jq- 2v j)} + R ( £ )- 

i<j ( q,l)€S Kql 

(17) 

To estimate the decay rate of the function V(x(t)), we use the following 
lemma. 


Lemma 3 Let V{x) be a function of class C 2 (D ) such that inequalities 

Pi\\ x \\ 2 <V(x) < (3 2 \\x\\ 2 , /?i > 0, (18) 

ai V(x) < ||VH(x )|| 2 < a 2 V(x), ai > 0, (19) 

( 20 ) 


dx 2 


< I 1 


hold with some constants a\, a 2 , Pi, /3 2 , /a in a convex domain D C M ? 
If x : [0, e] D is a function such that 

x(e) = x(0) — eW(i(0)) + r e , x(0) 7 ^ 0, 

with some r £ G W 1 , then 

V(x(e)) < V(x(0)) fl- ai e+^Lt + ^ 


( 21 ) 


2/?i m(o) 11 2 7ftiii(o)n 


( 22 ) 
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Proof. Let us denote x° = x(0), y = -£VL(i°) + r £ , and apply Taylor’s 
theorem for the function V(x° + y) with the Lagrange form of the remainder: 


V(x° + y) = P(x°) + 


i— 1 


dV(x ) 
dxi 


-Vi 


1 A 


% d 2 V(x ) 

+ 2 ^ dxtdxi ViVj 


i,j=1 


, (23) 




where 6 1 G (0,1). By applying the Cauchy-Schwartz inequality to expan¬ 
sion (1231) and exploiting assumptions (fl9l) . ( 120 T) . we get the following esti¬ 
mate: 

V(x(e)) < V{x°)-e (l - ||VP(r°)|| 2 +(l+£/i)||VP(a; 0 )|| ||r £ || + ^||r £ || 2 < 


< ^1 - a x £ + — V(x°) + (1 + £y)^/a 2 v(x°)\\r £ \\ + < 


< V{x°) 


1 — Ol\£ T 


a 2 £ 2 y 


+ 


V^ 2 (l + £fj)\\r e 


V\\re\\ 2 ) 
2V(x°)) 


(24) 


if V(x°) 7 ^ 0. Then the application of estimate (fl 8 |) to (1241) yields ineqial- 
ity ([ 22 ] 1 ) . □ 

By using Lemmas [ 2 ] and [3] for 7 r £ -solutions of system (|T|) corresponding to 
a partition 7 t £ = {sj}j >o and control u = u £ (t, x), we prove Theorem Q] 


4 Proof of the main result 

In order to prove Theorem [H let us first prove auxiliary lemmas. 

Proof of Lemma [lj. Let us enumerate the elements of S in ([2]) as 

1(^1; jl)i (*2)j2)) ■■■ j (^n —mi jn— m)} 


and introduce the n x n-matrix 

A ( x ) = ■■■Jm(x),[f il J jl ](x), ~;[fin- m J jn - m (x)]) ■ 

As the vector holds fi(x) satisfy the rank condition (J 2 j), there is a closed 
bounded domain O C D, 0 £ int D such that the map 

$(x) = -24 _1 (i)W(i) (25) 
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is defined for each x G fl and continuous. To study the solvability of equa¬ 
tions m, we introduce new variables 


“ 4t rhj 




and rewrite system (J5J) as follows 

F*x{£,) T G x {£) 0, £ (^l) v mi , 

where 

iE(£)=£-<f>(x), 

2 
-i 

£ 


-G*(0 = vpyA 
i,j =1 


- 1 / \ 
f x)-^ —-/*(a:)+ 


chc 








E 






A ( x )lfi,fj]( x ) 


i<j y ( q,i)£S V ( q,j)&S V ^ 

Here the integer constants k^ may be chosen as 

kiiji = I; ^*272 = 2, 


n—mjn—m 


= n — m. 


(26) 


If 5 = (t>i , v m , £a hjl , ...,£a in _ mJn _ m ) / is a solution of system (j2EJ) for given 
x G M” and £ > 0, then the components of a solution of equations (J5]) are 


vi,v 2 ,...,v m , dij = 2 y 7rfcjj | dij | sign , (i,j) G S', (27) 

with = kij if dij > 0 and k^ = —k^ otherwise. Tims, the solvability of 
system (JHJ) is reduced to the study of equation (1261) . 

Note that $(x) = 0 for x = 0 as V(x) is positive dehnite, so £ = 0 is a 
solution of equation (1261) for x — 0. 

To prove the existence of solutions for equation ([26]) . we find a A > 0 
and show that the degree of a continuous map 


£ G S p f—)■ 


ffr(£) + <?*(£) 

m)+G x m 


eSi, 


(28) 


is equal to 1, under a suitable choice of p > 0 depending on x if 0 < ||x|| < A, 
Ba( 0) C fh Here the spheres 


S p = {£ G M” I lien = P j and S, = {£ G R" | ||£|| = 1} 
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are oriented as (n — l)-spheres in M n . 

As is compact then there exist positive constants M 0 , Mi, and L such 
that 


\\A-\x)\\<M 0 , \\fi(x)\\ <Mi, 


dfi(x) 


dx 


< L, i — 1, 2, m, \/x G fh 


If ||£|| G S p then the Cauchy-Schwartz and triangle inequalities yield 


110.(011 <eM 

i<j 


\viVj\ + ^=-(n - m) 3/4 ||a|| 1/2 ^ (|ui| + \vj\) < 

i<j 


7T 


4 n 5 / 4 

where M = LM 0 Mi. We have also exploited formula (127]) together with the 
following properties of components of £ here: 



< £Mn||n|| 2 + 


4 eMn 9 ^ , 


7r 


I 1 / 2 < Mn\Jep 3 


Ml < IIS Nl<-ll£l|. 


(30) 


Then the maps _F X (£) + G x (£) and £ are homotopic on S p provided that 


ll£ll> 11^)11+ l|Gx(£)||, V£e5, (31) 

To satisfy condition fj3Tj) . we observe that ||W(a;)|| = 0(||x||) in a neigh¬ 
borhood of x = 0 for a positive definite function V(x). Hence, there exist 
positive constants A and A such that 


\Mx)W <il>\\x\l Vx G -Ba(O)- (32) 


By taking into account inequalities (1291) and (J32J), we conclude that condi¬ 
tion (j3TT) is satisfied for x G 0) if 


p > V’lkll + 


__ f4n 5 / 4 



The function 0(p) = p—Mny/ep 3 j = P+0(p 3 ^ 2 ) takes positive 

values for p e (0,p 0 ); where 


VPd = r 0 


27m 5 / 4 / f 7r \ 

VV 1 + 4Mn 7 / 2 “ ; 
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is the positive root of the following equation 


r 2 o + 


4 n 5 / 4 


£71 


r o - 


Mne 


0. 


Then we choose any r G (0, r 0 ) and check that 

0(p) > Kp for 0 < p < f 2 , 


where 

f 4?r 5 / 4 1 

A' = 1 — y/eMnf 3 | + fi/ij > 0. 

Let us take 

A f _ Kr 2 ) 

A = min < A, —— > > 0 

and observe that, for any x:0< ||a:|| <A, condition (l3l| holds if 

P=jf\\ x \\- ( 33 ) 

Homotopic equivalence of the maps F x (£) + G x (£) and ^ on S p , ensured 
by (j3Tj) . implies that the degree of the map fl28j) is equal to 1, i.e. to the degree 
of the map ||||j- : S p —> S\. By exploiting the degree principle (see, e.g., [IB]), 
we conclude that there exists a £ G B p ( 0) such that F x (£) + G X (£,) = 0, which 
means the existence of a solution to equation (j5j) according to formulas (j27l) 
if x G -Ba(O). Then estimates (TlTh follow from inequalities (j30l) and fj33jl . □ 
A useful a priori estimate of the solutions of system (JT|) is given by the 
following lemma. 

Lemma 4 Let x(t) GhC 0 < t < r, be a solution of system (P) with 
the control u G C[0,r], and let 

II fi(x') - fi(x ")|| < L\\x' -x"\\, W,x" G D, i = 1,2 

Then 

\\x(t)-x(0)\\<^-(e LUt -1), t G [0,r], (34) 

where 

m 

M = max ||/i(x(0))||, U — max \ui(t)\. 

Ki<m 0 <t<r L ' 

i=l 
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Proof. By differentiating the function w(t) = ||a;(f) — x(0)|| along the tra¬ 
jectory of system flTD, we get 


d 9/ . 

Jt w (t) = 2 


m \ 

x{t) -x( 0 ),^Ui(t)fi(x(t))J 


< 


< 2 Uw(t) max || fi(x(t)) - fi(x(0)) + /i(x(0))|| < 2 Uw(t)(Lw(t) + M), 

Ki<m 


SO, 


w(t) < U(Lw(t ) + M), t > 0. (35) 

We solve the comparison equation for differential inequality (l35jl to obtain 
the following estimate (cf. [14) Chap. Ill]): 


w(t) <^-(e LUt - 1), t G [0, r\. 


This proves estimate (134]) . □ 

Now we use Lemma ED to prove Lemma [2j 

Proof of Lemmo H For a solution x(t) of differential equation ([1]) with 
the initial condition x( 0 ) = x° and control u € C[ 0 , r], we represent the 
coordinates of A x(t) — x(t) — x° by the following integral equations: 


Ax i 


m r. 7 

w = £/ 

i=l 


Ui(t)f i 



Uj( s )fj(x° + A x(s))ds ) dt = 




2=1 


dfj(£(s)) 

dx 


Ax(s) ) ds + 


1 / d 2 f ik (r](t )) 

2 


0 2 Ax(t), Ax(f) j }> dt, k — 1,2,..., n. 


(36) 


Expression (1361) is obtained by Taylor’s theorem with the Lagrange form of 
the remainder, 0 < ||£(s) — x°|| < ||Aa;(s)||, 0 < || r)(t) — x°|| < ||Ax(t)||, for 
0 < s < t < t. Comparing formula (TT3|) with (136|) . we get 


Rk(r) = 


m 

E 

hi = 1 


dfik(x° 
dx 



Ui(t)Uj(s) 


0 Jo 


dfj^(s)) 

dx 


A x(s) ds dt+ 


+ 1E jT Al(( ), Mt) ) Mi(( ) dt. (37) 
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We use estimate (j5H) from Lemma Q] and the triangle inequality together 
with the Cauchy-Schwartz inequality to evaluate the Euclidean norm of the 
vector R(t) = (Ri(r ),..., R n (r)y in (l3?j) : 

\\R(t)\\ < L 2 U 2 j f ||Ax(s)|| dsdt+ f < 

Jo Jo 2 J o 

< f |e“ T - \ (( LUt + l) 2 + 1 )]+ H f I f' { C" T - 2 ) 2 + 2ic7 T - l} . 

(38) 

The right-hand side of formula ([131) is obtained as the Taylor expansion of 
formula (155]) with respect to Ur. □ 

Proof of Theorem [0 Let us denote D 0 = B Po (0) C D and choose a 
positive number e < e 0 such that all solutions x{t) of system (P) are well 
defined on t G [0, e] for each e G (0, i], provided that x° = x(0) G D 0 and the 
control u = u £ (t,x °) is given by formula (J3J) with parameters 

Vi = v-(x°), ciji = ctj^x 0 ), kji = k £ t fx°) 

obtained from algebraic system (J5]) . 

We define 

M= sup \\fi(x)\\,d= inf ||x — y\\ > 0. (39) 

xGDo, l<i<m x£D 0 ,y£OD 

If D = M n then we take d = +oo and i = £o, otherwise e < £ 0 is obtained as 
a positive solution of the inequality 

M (, e LU(i)£ _ 1) < (40) 

where L is given in condition (181). 

m 

U(e)= sup ^2\u £ (t,x)\, 
te[o,e],xeD 0 i=1 

and uf(t, x) are given by flU) . Condition (jUJ) implies that U{i)i —> 0 as £ —> 0. 
Thus, the set of solutions i G (0, e 0 ] °f inequality (140]) is not empty. Let e be 
such a solution, then from inequality ()40l) and Lemma [4] it follows that 

||x(f) — x°|| < d, t G [0, er] , (41) 
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for each solution x(t) of system ([Tj) with x° G D 0 and u = u £ (t, x °), e G (0, i]. 
Inequality (FXTj) means that x(t) G D for t G [0,e]. 

Let V (x) be a function that satishes conditions (J7|) . We introduce level 
sets 

L c = {x G D | V(x) < c} 


and define 


c 0 = inf V(x), p= inf ||x 

xGD\Dq x£D\L C0 


It is easy to see that Co and p < po are positive numbers as V (x) is positive 
definite. By the construction, 


B p ( 0) C L co C D 0 and L c C L co 


for each c < Co- 

The next step is to show that, if £ > 0 is small enough, then there exists 
a positive a — a(e) < 1 such that 

V{x{e)) < (l-(x)I/(x 0 ), (42) 

for any solution x(t) of system (JTj) with the initial data x° G L co and the 
control u = u e (t,x° ) given by (J3J). 

As V G C 2 (D) is positive def ini te then VV^(O) = 0, and Taylor’s theorem 
implies the following inequality: 

V(x) < f3 2 \\x\\ 2 , Vi G h 0 , (43) 


where 

d 2 V(x) 
dx 2 

is finite by Weierstrass’s theorem due to the compactness of D 0 . By applying 
similar argumentation to the function || VId(a;)|| 2 , we conclude that 

||VW(m)|| 2 < d 2 ||x|| 2 , ViGh 0l 

with some positive constant oti. Because of conditions (JT|) , it follows that 

\\VV(x)\\ 2 <a 2 V(x), Vx G D (h (44) 

where a 2 = d- 2 /> 0. Inequalities (|7J), (l43]h and (jTTl) imply that all con¬ 
ditions of Lemma [3] are satished in Do if x(t) (0 < t < e) is a solution of 
system (JTj) with the control u = u e (t, x°), x° G D 0 . 


2/? 2 = p = sup 

x£Dq 
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In order to satisfy condition (1421) . it suffices to assume that 

aiE -—~mw ——- • 6 oV{} (48) 

because of Lemma [3] Here the remainder R(e) of the Volterra series can be 
estimated by Lemma [2] as follows: 


\\R(e)\\ < HW 3 (x°)e 3 for x° G D 0 and W(x°)e < 1. 


(46) 


Here H is a positive constant, 


m 


W(x")= sup yj|<(i,i°)|, (47) 

te[0,er] i=1 

and u £ (t,x °) is given by (JH). Condition Q together with representation (]4|) 
implies that 

W{x°) < C||a: 0 || 1/3 £- 2/3 (48) 

with some positive constant C for all x Q G D 0 . 

Estimates (1461) and (148]) imply that condition (145]) holds if 

a 2 fie hH 2 C 6 £ 2 y^(l + Efi)HC 3 e 


(X\ — 


2/3, 


Vfii 


> a. 


(49) 


Here a = o/e is a positive number. As is positive, we conclude that there 
exist E m ax > 0 and <7 > 0 such that inequality (j49j) holds for all £ G (0, e ma J. 
Without loss of generality we suppose that such Emax corresponds to the 
assumption of formula (]46]h i.e. W(x°)£ m ax < 1 for all x° G D 0 . Thus 
we have proved that condition (|45|) is satisfied for each e G (0,£ max \ with 
a = <t(e) = min((je, 1). Let us define e = min(e, £ max ), where i is a positive 
solution of inequality (HOjl . Then inequality (l42|) holds for any £ G (0, e] with 
a = cr(£) < 1 provided that x° G L c o. 

If x° G B p { 0), £ G (0, e] and u £ (t,x) is given by formula (j3J), then the 
corresponding 7T e -solution of system (jT]) x(t) is well-defined: 

x(ne) G L co C D 0 for n = 0,1, 2,..., 

and x(t) G D for all t > 0 because of inequality (FXTf . By iterating inequal¬ 
ity (1^2]) for x° G B p ( 0) C L Co , we conclude that 



for t — 0, £, 2£,..., 


(50) 
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where 


A 


ln(l — a) 

Ye 


> 0 


if cr < 1, 


and A is an arbitrary positive number if a = 1. For an arbitrary t > 0, we 
denote the integer part of 7 as N — [j] and denote t — t — Ne < e. Then 
we apply inequality (15U|) together with Lemma 3] to estimate x(t): 


|x(i)|| = || x(t) — x(Ne) + o;(iVer)|| < ||x(lV£)|| + ||x(t) — x(iVe)|| < 


< * —\\x 0 \\e~ XNe + — (e LW ^ N£ » £ - 1) , 


(51) 


where L, M, and W{x) are defined in (JHJ), ([39]) . and (147)) . respectively. Esti¬ 
mates (148|) and (150|) imply the following asymptotic representation: 


W(x(Ne)) = 0(||a;(iV£)|| 1 / 3 ) = 0(e ~ XN£/3 ) as N -> + 00 . 


Then if follows from inequality ()5T]) that 

||x(£)|| = 0(e~ xt ) as f-> +00 

with A = A/3>0. By using formulas flTTl) and (l48]h we conclude that 
\\u £ (t,x(t))\\ = 0(e~ xt ^ 3 ) as t —> + 00 . 


□ 

5 Stabilization of the Brockett integrator 

Consider the control system known as the Brockett integrator |9j: 


Xi = Ui, X 2 = U 2 , x 3 = U 1 X 2 - U 2 x 1 , (52) 

where x — (aq, x 2 , x 3 )' G R 3 is the state and u = (ui,u 2 )' G R 2 is the control. 
The stabilization problem for system (1521) has been has been the subject of 
many publications over the past three decades (see, e.g., the book [5] and 
references therein). I 11 particular, it is shown that system (1521) can be ex¬ 
ponentially stabilized by a time-invariant feedback law for the initial values 
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in some open and dense set fl 7 ^ R 3 , 0 ^ int Q |H] - In this section, we con¬ 
struct a time-varying feedback law explicitly in order to stabilize system (1521) 
exponentially for all initial data. 

System (152|) satisfies the rank condition of form (J2J) with S = {( 1 , 2 )}, 
span{ /i (x), f 2 (x), [fi, f 2 ]ix)} = R 3 for each x G R 3 , 


where the vector fields are fx{x) = ( 1 , 0 , a? 2 ) / , fi{x) = ( 0 , 1 , —Xi)', 

0 . 0 ,- 2 )'. 

The family of controls (fl 6 |) takes the form 


ui(t) = Vi + a i2 cos {^^t) , 

u 2 (t) = v 2 + a i2 sin , fc 12 G Z \ {0}. 


(53) 


For an arbitrary initial value a : 0 = (a;) 1 , a; 2 , x°)' G R 3 at t = 0, the solution 
x(t) of system ([52]) with controls (T53l) is represented by (TlTl) as follows: 


xi(e) = a:° + ev\, x 2 {e) = a;° + £v 2 , 

E 2 

x 3 (e) =xl + e (vix% - v 2 x°fi - —-—a 12 (a 12 - 2vf). (54) 

Z7T/Ci2 

Note that representation (1541) is exact (i.e. the higher order terms R(e) in 
the Volterra expansion vanish) as system (|52l) is nilpotent. This implies the 
following lemma. 


Lemma 5 For arbitrary x° = (x^, x®, x®) 1 G R 3 , a : 1 = (x{, a; 2 , a^)' G R 3 , and 
£ > 0, define the controls Ui = ufit) and u 2 = u 2 (t) by formulas (153]) with 


v v = 


(V* 1 _ /y»0 

Jb 1 Jb 


v 2 = 



£ 


a 12 = ——— ± - \J(x\ — X ?) 2 + 2iTki 2 (xl — X 3 + x\x^ — 

Then the corresponding solution x(t) of system (1521) with initial data x(0) = 
x° satisfies the condition x(e) = x 1 . 
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To solve the stabilization problem for system (15211 . consider a Lyapunov 
function candidate 

= \{ x l + A + xl). 

Following the approach of Theorem [U we define a time-varying feedback 
control to approximate the gradient flow of —'VV(x) by trajectories of sys¬ 
tem (1521) : 

/ 2irk(x ) 


u\(t, x) = r?i(T) + a(x) cos 


V £ 


-t 


where 


. . . . . . ( 2nk(x) 

2 (t, x) = v 2 {x) + a{x) sm I - 1 


v\(x) = —xi, v 2 (x) = —x 2 , k(x) = signx 3 , 

a(x) = ( ~ Xl ± + ~f~ 0, 

0, x 3 = 0. 


(55) 

(56) 


Without loss of generality, we may assume any integer value for k{x) if x 3 = 0. 

By Theorem [TJ the feedback control (l55l) - (l56l) ensures exponential stabil¬ 
ity of the equilibrium x — 0 in the sense of 7r £ -solutions, provided that £ > 0 
is small enough. 


6 Simulation results 

In this section, we perform numerical integration of the closed-loop sys¬ 
tem (l52l) with the feedback law u = u(t,x(t)) of form (l55l) - (j56l) . Trajectories 
of this system are shown in Fig. [Hand [2] for £ = 1 and the following initial 
conditions: 

®i(0) = x 2 (0) = 0, x 3 (0) = 1 (Fig. □), 

®i(0) = x 2 (0) = x 3 (0) = 1 (Fig. ED . 

These simulation results show that the feedback law ([55jl -(l56 jl steers the 
Brockett integrator to the origin not only in the sense of 7T e -solutions (as 
stated in Theorem [1]), but also in the sense of classical solutions. 
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Figure 1: Trajectory of the closed-loop system f[52jl . (l55jl . (I56j) for Xi(0) 
x 2 (0) = 0, x 3 (0) = 1. 



Figure 2: Trajectory of the closed-loop system (l52]h (l55]h (I56j) for Xi(0) 
^2(0) = x 3 (0) = 1. 
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7 Conclusion 


In this paper, a family of time-dependent trigonometric polynomials with 
coefficients depending on the state has been constructed to stabilize the 
equilibrium of a nonholonomic system. These coefficients are obtained by 
solving an auxiliary system of quadratic algebraic equations involving the 
gradient of a Lyapunov function. An important feature of this work relics on 
the proof of the solvability of such a system for an arbitrary dimension of the 
state space provided that the Lie algebra rank condition is satisfied with first 
order Lie brackets. It should be emphasized that this result is heavily based 
on the degree principle as the implicit function theorem is not applicable for 
a non-differentiable function G x (£) in Lemma [TJ. 

Another important remark is that our design scheme produces small con¬ 
trols u £ (t,x) for small values of ||x||, and the frequencies of the sine and 
cosine functions are constant for each fixed e > 0. This feature differs from 
the approach to the motion planning problem that uses a sequence of high- 
amplitude highly oscillating open-loop controls (see m inn- 

The proof of Theorem |T] is considered as an extension of Lyapunov’s direct 
method, where the decay condition for a Lyapunov function is guaranteed by 
exploiting the Volterra expansion instead of using the time derivative along 
the trajectories. Although the exponential stability result is established for 
7r £ -solutions under a sampling strategy, simulation results demonstrate the 
convergence of classical solutions of the closed-loop system to its equilibrium. 
Thus, the question of the limit behavior of classical (or Caratheodory) so¬ 
lutions of system (JT]) with the feedback control (J4]) remains open for further 
theoretical studies. 
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